The Schwinger mechanism revisited 
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The vacuum persistence probability, P vac (t), for a system of charged fermions in a fixed, external, 
and spatially homogeneous electric field, was derived long ago by Schwinger; w = — lag[P vac (t)]/Vt 
has often been identified as the rate at which fermion-antifermion pairs are produced per unit volume 
due to the electric field. In this paper, we separately compute exact expressions for both w and for 
the rate of fermion-antifermion pair production per unit volume, V, and show that they differ. While 
w is given by the standard Schwinger mechanism result of w = ^§3- X/nLi ex P ) 1 ^ ne 

pair production rate, F, is just the first term of that series. Our calculation is done for a system 
with periodic boundary conditions in the Ao = gauge but the result should hold for any consistent 
set of boundary conditions. We discuss, the physical reason why the rates w and F differ. 

PACS numbers: 



I. INTRODUCTION 

Fermion-antifermion pair production from a static clas- 
sical electric field, known as the Schwinger mechanism, 
has given rise to a vast literature, since its formulation 
in 1951 [ij. Invoked to gain insights on topics as diverse 
as the string breaking rate in and on black 

hole physics |5j, this mechanism has become a textbook 
topic in quantum field theory 0. Topics such as back 
reaction [loj and finite size effects^ have been addressed. 
In a classic paper, Schwinger exactly calculated the rate 
at which the vacuum decays due to pair production in the 
external field. If the electric field is treated classically — 
i.e. the formal limit of q — > 0, E — > 00 with qE fixed — the 
vacuum persistence probability is[l[: 



P vac (t) ee |(vac|C/(t)|vac)| 2 =cxp(- W yt) 



with w 
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E; 

n=l 



exp 



qE 



(1) 
(2) 



where V is the spatial volume of the system and w is the 
rate of vacuum decay per unit volume. 

While the Schwinger formula of Eq. ([2]) is very well 
known, the Schwinger mechanism has often not been 
fully appreciated in an essential way. The quantity w 
in Eq. ^ appears to have a very natural interpreta- 
tion as the rate of production of pairs per unit volume. 
Schwinger suggested this interpretation in his original pa- 
per. There is a plausible, if heuristic, argument in its 
favorJ^Q. Start by considering a more general situation 
in which the pair production rate can vary in space and 
time. The vacuum persistence probability can then be 
written as 



-Jd 4 



'j(x) 



(3) 



Next approximate the integral as a discrete sum over 
space-time cells of volume Avi points centered at points 
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Xi- The vacuum persistence probability then becomes 
the limit of 



-Pvac = lim TT 1 



J{1 - w{xi)Avi) . (4) 



Itzykson and ZuberQ note that this is the form expected 
from an independent contribution to the vacuum persis- 
tence coming from each cell if w(xi) is the local rate of 
pair production. This is taken to confirm Schwinger's 
interpretation of w. This argument is plausible; it is 
not surprising that a considerable body of literature has 
adopted Schwinger's interpretation and used w to com- 
pute particle production rates. 

However, as it happens, the interpretation is not cor- 
rect: in general w does not give the rate of pair creation. 
As a logical matter the rate characterizing the vacuum's 
decay (by the production of the first charged pair), is not 
necessarily the same as the continuous rate of pair pro- 
duction. While the argument given above is highly plau- 
sible, it is also heuristic. A more compelling approach 
is via a comparison of a direct computation of the pair 
production rate with w, the rate in Eq. @. Nikishov di- 
rectly computed the pair production rate per unit volume 
long agoQ and found it to be 
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Remarkably it does not agree with w: the entire rate is 
given by the first term in the series for w. 

Since w is still commonly confused with the pair pro- 
duction rate, it is useful to rederive Nikishov's result in 
a physically transparent context which clearly illustrates 
why r rather than w is the rate of pair creation per unit 
volume. This will first be done in the context of (1+1) 
dimensions where the issues are particularly clear; the 
generalization to (3+1) dimensions is straightforward. 

Before embarking on the calculation, it is useful to dis- 
cuss briefly why one might expect T (the pair production 
rate per unit volume) to differ from w (the rate char- 
acterizing the rate of vacuum decay per unit vacuum). 
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We begin with the trivial observation that pairs are cre- 
ated in distinct modes. If the probability for any given 
mode to create a pair is constant in time (and thus un- 
corrected with each other) , one expects that total rate of 
pair creation to equal w. However, if there are temporal 
correlations between likelihood for decay at various times 
for the various modes, then this need not be the case. 

The precise time that a pair is created is ambiguous; 
however, once a pair is created and the fermion and an- 
tifermion become well separated, they each accelerate in 
the electric field. Thus particles with large (mechani- 
cal) momentum are likely to have been created earlier 
than particles with small momentum. To the extent that 
the natural way to characterize modes is in momentum 
space, there are very large temporal correlations in pair 
creation. As will be seen in detail in the explicit cal- 
culations, momentum space is the natural way to char- 
acterize modes for this problem and it is natural that 
r differs from w. This makes clear why the argument 
of Ref. that w is the pair creation rate breaks down: 
there it is assumed that the rate of pair production at dis- 
tinct space-time points are independent. Implicit in this 
is the notion that the pairs are created at well-defined 
points in space-time. However, if the pairs are created at 
well-defined (or nearly well-defined) momenta they are 
delocalized in position. 



II. PAIR CREATION IN (1+1) DIMENSIONS 

Before embarking on the original problem, it is use- 
ful to consider the analogous problem in (1+1) dimen- 
sions. The (1+1) dimensional case is somewhat sim- 
pler and the issues are more straightforward. Moreover, 
the (1+1) dimension result serves as an important input 
into the (3+1) dimensional calculation. One can adapt 
Schwinger's calculation [l| for smaller dimensions in 
(1+1) dimensions P vac {t) is 

P^(t) = exp(-u, 1 + 1 Lt) 



with w 1 




where L is the length of the system. 

The imposition of appropriate boundary conditions in 
both space and time is crucial in the treatment of this 
problem. Accordingly, it is critical to compute both w 
and r, with the same boundary conditions so they can 
be directly compared with each-other. We adopt the fol- 
lowing strategy in choosing boundary conditions in space 
and time. 

The most natural way to identify the number of pairs 
created is to formulate the problem in terms of an elec- 
tric field which is switched on for a fixed period and then 
switched off. One begins with the system in its vacuum 



state and the electric field turned off. The electric field 
is turned on at t = and left on until t = T and then 
switched off. With the field now off, one can unambigu- 
ously use a free particle basis to count the number of 
pairs. The act of turning the system on and off can yield 
transient effects; thus one needs to consider the limit of 
large T (T > m _1 , T > qE^ 1 / 2 ) so that the steady 
state Schwingcr mechanism dominates. 

To simplify the computation it is useful to keep the 
electric field a constant over the entire system while mak- 
ing the system of finite length. In this case the modes 
become discrete and are easily counted. The large L limit 
can be taken at the end. Continuous acceleration of par- 
ticles over extended time plays a key role in the analysis. 
Thus it is important to choose spatial boundary condi- 
tions which allow a charged particle to continue accelerat- 
ing as it reaches the end of the system. Dirichlct bound- 
ary conditions are thus inappropriate as a particle strik- 
ing the end of the system will bounce back and, instead 
of continuing to accelerate, will then decelerate. Contin- 
uous acceleration over long times can be achieved simply 
via the implementation of periodic boundary conditions. 
However, periodic boundary conditions are sensible only 
when the Hamiltonian is invariant under translations by 
L. Because the gauge choice Aq = ensures transla- 
tional invariance, we work in this gauge. Of course, our 
results must be gauge invariant and we can work in any 
gauge. However, if we choose to adopt a different gauge 
the boundary conditions needed to ensure continuous ac- 
celeration become time dependent. Ultimately we can 
check whether this choice of boundary condition (and its 
associated gauge choice) is sensible by comparing the de- 
rived vacuum persistence probability with the Schwingcr 
result of Eq. © . The Dirac equation becomes 

(a(-id x - qA x (t)) + f3m)ijj(x,t) = id t ip(x,t) 

{0 for t < 

—Et for < t < T 
-ET for T < t 

ip{x,t) = ip(x + L,t) (7) 

where ip is a two-dimensional spinor and a and j3 are 
the appropriate two-dimensional Dirac matrices. There 
exists a complete set of a solutions of the form 

^(x,t) = e^ x X k(t) 

Pk = k — (integer k) 

h k (t) = a(pk - qA x (t)) + (3m 
hk{t)Xk{t)=id tX k{t) (8) 

where the restriction to integer k comes from the peri- 
odic boundary conditions. Finally, it is convenient to 
introduce a unitary transformation, Uk(t), into a basis in 
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which h k is diagonal 

Xk(t) 
e k (t) 

Ut(t)h k (t)U k {t) 



Pk - 

<*(*) 




q A{t)f 


-c fc (t) 



For convenience we choose T such that T 



JT 



(9) 



with 



r = and j a positive integer. It is easy to sec that 
e k (T) = e k+ j(0): the spectrum of the system after E is 
turned off is identical to the original spectrum, although 
the mode labels have changed. This is hardly surprising, 
as the system for t > T is simply that of a free par- 
ticle in the absence of an electric field — exactly as it is 
for t < 0. The restriction to integer j ensures that the 
boundary conditions are the same. At t — the system 
is in its vacuum state with all of its negative energy levels 
filled. Thus the appropriate boundary condition for the 
fermionic modes is $jT(0) = (0, 1). 

The equation of motion for = (<P k {t), <t> k (t)) is 



obtained from Eqs. 

by 

I qEm 



and © and for < t < T is given 




id* 



Kit) 



(10) 



By construction, \4> k (T)\ 2 is the probability for the posi- 
tive energy state of the fc th mode to be occupied after the 
field has been switched off; it represents the probability 
of pair creation for this mode. The expectation value for 
the total number of pairs produced, TV, and the vacuum 
persistence probability arc thus given by 



(^cn>=£i<^cr)i 2 

k 

P^ 1 (T) = il(l-\ti(T)\ 



111) 



Zenerfli"l|; exp ^— is simply the Landau-Zcncr tran- 
sition probability which is valid when the initial and fi- 
nal levels are well separated on the scale of m. Pro- 
vided that qET 3> to the correction to Eq. (TT3"|) is 
small except for a small number of modes. Inserting 
this form into Eqs. (fTTj) and (fTS]) and using the identity 
log(l - 6(y)6(z)x) = 9(y)6(z) log(l - x) yields 



(N(T)) = 0(-Pk)6(Pk + qET)e~ 



x [l + O 



qET 



(14) 



log (P4+ X (T)) = J2 0{-Pk)0{Pk + QET) log (l - e 



x l + O 



qET 



(15) 



where the scale of the correction term reflects the fraction 
of modes for which the correction to the leading term in 
Eq. p3[) is non-negligible. From the definition of pk in 
Eq. © one sees that J2 k 6{-p k )6(p k + qET) = 2f^. 
Using this fact, and taking the long time limit to remove 
transient effects, we see that P^+ 1 (t) in Eq. (fTS]) exactly 
reproduces the standard result for the vacuum persis- 
tence probability in Eq. ([6]), indicating that our boundary 
conditions were chosen consistently. 

Further, the rate of pairs produced by the electric field 
in the long time limit is, 



-l+i 



(iV) 
LT 



qE 
2tt 



exp 



qE 



(16) 



Critically, we see that the rate of pairs produced 

per unit length per unit time, is not u> 1+1 but rather the 
first term in the series as expected from Ref. Q- 



The massless limit 



exp Vlog(l -|</>+ (T)| 



(12) 



To proceed further we need solutions of Eq. (|10[) . An 
exact solution can be obtained in terms of parabolic 
cylinder functions with complex arguments. The form 
is cumbersome and will not be given here. The impor- 
tant point here is the asymptotic behavior valid when 
\p k \ 3> m, and \p k + qET\ 3> m, in which case the solu- 
tion reduces to 



I^COI 2 = 0{-p k )6{ Pk + qET) exp 



l + O 



qE 



\p k y\p k +qET\ 



(13) 



The form of Eq. (|13p is not surprising. The equation 
of motion for the two level system in Eq. ([S]) is pre- 
cisely of the form considered long ago by Landau and 



The m — > limit in (1+1) dimensions illuminates the 
essential issues quite clearly. From Eq. © it is clear that 
diverges as to — * 0. Thus, the vacuum persistence 
probability becomes zero. If the interpretation of w as 
the pair production rate were correct this would mean 
that the rate of pair production per unit length would, 
perversely, be infinite for massless particles. However, 
the particle production rate per unit length is given by 
rl^lin = and is finite. 

It is easy to visualize this physically. With the bound- 
ary conditions used here, the energies of the modes for 
to = are given by +{p k + qEt) as shown in Fig. [1] 
Moreover the occupation numbers are preserved by the 
equation of motion; only the mass term induces transi- 
tions. There is a subtlety when p k + qEt = 0; at which 
point the positive and negative energies cross. Since <j) + 
indicates the amplitude for the positive energy solution, 
the labels (j> + and cf>~ switch at p k + qEt = 0. (Note 
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1 2 3 4 5 




1 2 3 4 5 



t/T 

FIG. 1: Energy levels for massless (1+1) charged fermions 
in a constant electric field as a function of time in units of 
r = . The dotted lines correspond to empty-initially pos- 
itive energy-levels and the solid lines to filled-initially nega- 
tive energy-levels; the dashed lines are the half-filled levels of 
the zero-mode. The solid circles correspond to particles and 
the open circles to antiparticles; the half-filled circles corre- 
spond to the creation of a particle (antiparticle) with 50% 
probability 

that when any finite mass is put in, the level crossings 
are avoided. However, as m — > the Landau-Zener prob- 
ability which the occupation number switches from the 
negative to positive level goes to unity and one repro- 
duces the zero mass result.) 

Initially, the system is in the vacuum state with all 
negative energy levels filled. Thus except for the special 
case of k = 0, |^ + (i)| 2 is given by 9(-p k )0(qEt +Pk), 
exactly as one expects from Eq. (|T3|) . The k = modes 
are special in that they have exactly zero energy at t = 0; 
we take them to be half occupied. Consider the occupa- 
tion of positive energy levels as time increases. As seen in 
Fig. [1] when t increases by r = exactly one new level 
crosses from negative to positive and this corresponds to 
the creation of a pair. At t = jr one has created on aver- 
age j — 1/2 pairs (the 1/2 coming form the 1/2 filled k = 
mode): the number of pairs increases linearly with time 
at a rate of 1/r = a j—, precisely the rate in Eq. (fT6)) . 
However, for j > 2 the vacuum persistence probability 
is strictly zero; at j = 2 we know with unit probabil- 
ity that a pair has been created in the k = — 1 mode. 
This cleanly illustrates both the distinction between T 
and w and the critical role played by temporal correla- 
tions between modes. It is the fact that pairs are created 
in different modes at different times which accounts for 
the difference. 

A natural way to visualize the massless case in (1+1) 
dimensions is to think of the (1+1) dimensional "uni- 
verse" as a circular loop of radius R = L/(2n). Suppose 
further that this loop is imbedded in a (3+1) dimensional 
world with a time- varying magnetic flux perpendicular to 
the loop and threading through its center. This flux is 
turned on at t — 0, after which it linearly increases in 



time, for a time t = jr after which it is held fixed; r is 
defined as the amount of time for one flux quantum to 
thread through. Because ^ = -Vx£, this is identical 
to an E field pointing tangent to the loop at all points, 
for a total time At = jr. The ground state of this system 
at t = is simply a filled Dirac sea of massless, charged 
particles — in just the same manner as the ground state 
of a normal wire represents a filled Fermi sea of massive, 
charged, particles (electrons). This system clearly has pe- 
riodic boundary conditions, and thus gives another phys- 
ical picture for understanding pair creation of charged 
particles in (1+1) dimensions. While the E— field is on, 
magnetic flux-quanta thread through the ring at a rate 
of 1/t. Since the system is circular, the natural variable 
labeling modes is the angular momentum, Lk (where Lk 
is related to the momentum in the previous formulation 
via Lk = Rpk)- After j flux quanta have thread through 
the loop, then Lk — > Lk+j ■ This both drives j empty lev- 
els into the Dirac sea and j filled levels out of the Dirac 
sea, creating j fermion-antifermion pairs. The rate is 1/r 
and is precisely that given in Eq. (|16p . 



B. A card game 

The distinction between the rate per length at which 
pairs over time, and the rate per length of vacuum 

decay, w 1+1 , is essentially one of classical probabilities; 
for this (l+l)-dimcnsional system this can be simply il- 
lustrated in the context of a card game. 

One might imagine that a casino has introduced a new 
game called "vacuum." The rules of the game are sim- 
ple: A very large number of decks of cards are shuffled 
together. One class of card is specified — this class could 
be rather general (e.g., all black cards or all diamonds), 
somewhat more specific (e.g., all sevens), or quite spe- 
cific (e.g., aces of spades). The pile of cards is said to 
be in the vacuum. Cards are then turned over at a fixed 
rate of 1/r. The pile remains in the vacuum until the 
first card of the specified class is turned over. Wagers 
are placed on how long the pile is in the vacuum and on 
how many cards of the specified class are turned over in 
a specified time. 

It is easy to see that if one continues to turn over cards 
at the rate of 1/r, the average number of cards in our 
specified class which are turned over after t = jr (where 
j is a positive integer) is (iV c i ass ) = ft/r where / is the 
fraction of the cards in the class (i.e., f = 1/2 if the 
class is black cards; / = 1/52 if the class consists of aces 
of spades). Thus, 7, the average rate of production for 
cards in the class, is given by 

7=-/- (17) 

T 

The probability of the pile remaining in the vacuum after 
t = jr is obviously given by P vac = (1 - /) J = (1 - /)'/ T . 
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This can be rewritten as 



P vac = exp(-wi) 

1 1 f n 

u =—iog(i-f) = -j2 — 

t t — ^ n 



(18) 



It is obvious that the structure of Eqs. p7|) and (|18|) are 
analogous to Eqs. (HHJ) and © with / playing the role of 



In the context of the game it is obvious why 7 and tu 
differ. Suppose as an extreme case one considers as the 
class all cards. The rate at which cards in the class are 
produced is obviously 1/r; after the first card has been 
turned over the probability that the pile is in the vacuum 
is clearly zero. This is the analog of the m — > limit. 



III. PAIR CREATION IN (3+1) DIMENSIONS 

Having illustrated the essential distinction between the 
vacuum decay rate, w, and the rate of particle produc- 
tion, r, in the (1+1) dimensional system, we turn to 
the case of (3+1) dimensions. For simplicity we spec- 
ify E = \E\z, and require solutions of the Dirac equa- 
tion to satisfy periodic boundary conditions along the 
transverse direction (making the spectrum discrete and 
counting explicit). In (3+1) dimensions there are two 
spin states for the fermions. However, in the energy 
basis analogous to Eq. the system is diagonal in 
spin; thus its sole effect is an overall factor of 2 in w 
and r. In direct analogy to Eq. (TlT)|) , the equation of 
motion for the mode specified by the spin state, s, and 
k = (k x , ky, k z ) = {n x — , n y — , n z — ) is 



qE 771 




(19) 



with 



►(*) 



yj(k z + qA z (t)) 2 + (j>l +p 2 y + m 2 ) 



the 



instantaneous energy of the mode k at time t. 
As in (1+1) dimensions, the initial conditions are 
(</>+ (0),</>+ (0)) = (0,1). Note that since (in this 

gauge) does not change the transverse momentum, it acts 
exactly like an additional mass-term in the Landau-Zcncr 
tunneling rate. Thus, for convenience, we will define 
mf, = m 2 + k\. 

All of the previous arguments leading up to Eq. (TT3"]) 



derive 

(N(T)) 
VT 



0(-k z )9(k z + qET) w 

2 2^ ~ ~ ~ ~ — e 



kx ,k z 

>qE L Z T 
2tt L 



L x L y L z T 



Evaluation of this double integral exactly reproduces the 
rate of pair production given in Eq. ^ . A similar calcu- 
lation of log[P vac (T)] ee -wVT yields, 



\og[P vac {T)} ee -wVT, 

= 2 6{-k z )0{k z + qET) \og[l-e 

qEL z T L x L 



7T 2tt 2tt 

\2 



1 dk x dk y log[l — e " E 



e " E 



(21) 



go through, after the substitution m 



ii. Thus, we 



Evaluation of the double integrals, and identification of 
this quantity with —wVT, yields the expression for w in 
Eq. ©• Exponentiating this rate, we derive the standard 
vacuum persistence probability in Eq. (TTJ). 



IV. DISCUSSION 

In summary, we have explicitly shown that the rate as- 
sociated with vacuum decay, w, and the rate of pair cre- 
ation, r, differ: T is given by the first term in the series 
for w. This result is in a very real sense quite well known: 
it was first derived by Nikishov Q quite long ago. It can 
be derived in various other elegant formulations [l2l [l3|. 
However, it is not as widely appreciated in the commu- 
nity as it should be: much of the literature in the field 
is still based on Schwinger's initial assumption that w 
gives the rate of pair production per unit volume. It is 
hoped that the physically transparent way that this pa- 
per shows that T differs from w will help clarify the issue. 
The massless case in (1+1) dimensions is particularly il- 
luminating. It is clear from Fig. [T] why the pair produc- 
tion rate remains finite even while the vacuum persistence 
probability goes to zero at finite times indicating an infi- 
nite value for w. While this calculation was with periodic 
boundary conditions in the Aq = gauge, we expect the 
result to hold generically for any consistent set of bound- 
ary conditions (i.e., those that reproduce the Schwingcr 
result for w). This will be investigated in future work. 

In practice, the distinction between if and T is expo- 
nentially small for weak fields. However, for strong fields, 
w exceeds I by a factor as large as ((2) — ^- ~ 1.64. 
More generally, in d space-time dimensions w d will ex- 
ceed T d by a factor of ((d/2) in the strong field limit. 
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